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ABSTRACT 

The barycenter of a massive black hole binary will lie outside the event horizon of the 
primary black hole for modest values of mass ratio and binary separation. Analagous 
to radial velocity shifts in stellar emission lines caused by the tug of planets, the radial 
velocity of the primary black hole around the barycenter can leave a tell-tale oscillation 
in the broad component of Fe Ka emission from accreting gas. Near-future X-ray tele¬ 
scopes such as Astro-H and Athena will have the energy resolution {SE/E < 10“^) to 
search nearby active galactic nuclei (AGN) for the presence of binaries with mass ratios 
q > 0.01, separated by several hundred gravitational radii. The general-relativistic and 
Lense-Thirring precession of the periapse of the secondary orbit imprints a detectable 
modulation on the oscillations. The lowest mass binaries in AGN will oscillate many 
times within typical X-ray exposures, leading to a broadening of the line wings and an 
over-estimate of black hole spin in these sources. Detection of periodic oscillations in 
the AGN line centroid energy will reveal a massive black hole binary close to merger 
and will provide an early warning of gravitational radiation emission. 

Key words: galaxies: active - accretion, accretion discs - black hole physics- rela¬ 
tivistic processes - techniques: radial velocities 


1 INTRODUCTION 


Galactic nuclei host su permassive black holes (> 10 ®Mq) 
llKormendv fc h 3 l2013l l and should be the sites of merg¬ 
ers of massive black holes (MBH). Binary black holes are 
expected to arise due to merging supermassive black holes 
(SMBH) throughout cosmic time in ACDM cosmology (e.g. 
ISpringel et al.l 20051 '). as well as due to intermediate mass 
black hole ( IMBH) formation in d isks around SMBH (e.g. 
Ibevinl [ 2 OO 7 I : iMcKernan et al.ll2012l '). Possible observational 
sign atures of these merging MBH are important to find 
fe.g. iMil osavlieyic &: Phinnevll200 ^: iBoedanoyic et al.ll2009l: 


Bodeet_^ 20101: iMcKernan et al.ll2oi^ 2014 : Roedig et al.l 

2 OI 4 I : iFarris et akl l2015l l as the search for binaries probes 


ever s maller separations (Comerford et al.ll2013l : iFord et aP 
I2OI4I : [Graham et al.l 12015 ). Here we discuss the possibility 
of detecting those massive black hole binaries in the local 
Universe closest to merger, via a radial velocity imprint in 
the broad FeKa line. We demonstrate the effect and its 
magnitude and we show that near future X-ray telescopes 
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with very fine energy resolution {Astro-H and Athena with 
SE/E < 10~® in the FeK band) will be able to detect MBH 
binaries near merger at mass ratios q > 0.01 and separations 
of several hundred gravitational radii. 


2 BINARY BLACK HOLE DYNAMICS 

A binary black hole will orbit its center of mass (barycenter), 
thereby adding a radial velocity component to line-emitting 
material immediately surrounding the primary SMBH. This 
is analagous to radial velocity shifts of stellar spectral lines 
as t he star orbits its barycent er due to planetary tugs 
fe.g. lOhta. Taruva fc Sutoll2005l 'l. By analogy, we can write 
the radial velocity compo nent (vrad) of the primary S MBH 
around the barycenter as llOhta. Taruva fc Sutdl2005l ') 

Vrad = Y sin(i)[sin(f -|- ZJ) -|- e sin(Zu)] (1) 

where q = m 2 /Mi is the mass ratio of the secondary BH 
(m 2 ) to the primary BH {Mi), Ub is the binary semi-major 
axis, = 1 — where e is the orbital eccentricity of the 
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secondary, ly = 2 tvITo^\j is the secondary mean motion or or¬ 
bital frequency, i is the inclination angle of the binary to the 
observers’ sight line {i = 90° is edge-on and yields maximal 
radial velocity components), / is the true anomaly or the 
angle made by the secondary in its orbit as measured from 
a line joining the focus and periapse, and a; is the argument 
of per iapsisu (see e.g. Figs 1 & 2 in lOhta. Taruva fc Sutol 
ll2005h for sketches of these orbital parameters). The inner 
accretion disk is bound to the primary so spectral lines orig¬ 
inating in the inner disk will oscillate around their rest cen¬ 
troid energies by Urad/c on the timescale of the secondary’s 
orbit (IYu fc Lull200ll : iBoedanovic et al]l2009l : ISesana et al.l 
I2OI2II . As binary separation shrinks, the disk can become 
strongly perturbed. In the limit of at of a few tens of Vg, 
we should expect the line profile to diffe r from the simple 
model below (see e.g. ISesana et al] (120121') : [McKernan et al.l 
(l2013tl and references in the latter for a discussion of the 
wide range of possible configurations including secondary 
disks, gaps and cavities). However most FeKa emission may 
originate within ~ 20 — 30rg of the primary and very tight 
binaries (at of tens of Vg) will be much rarer than those 
at wider separations (at of a few hundred Vg). Even if the 
line profile is perturbed, the periodic wobble will occur in 
the manner we outline below and may still be detectable. 
Eqn.[T] can be rewritten as a energy shift {SE) in a spectral 
line centroid energy {E) as 





Q sin(i). . 

TTq ~£ + ‘^) + ® sm(a;)] 


( 2 ) 


where rg = G{Mi + m' 2 )l<? is the gravitational radius 
of the binary and at is in units of Vg. Fig. [T] shows the 
maximum SEjE as a function of binary separation and 
mass ratio from eqn. [2] assuming circular orbits and edge- 
on inclination. Red horizontal dashed lines in Fig. [T] denote 
planned micro-calorimeter en ergy resolution for Astro-H 
fupper: [Takahashi et al.l (l201(ll 'l') and the micro-calorimeter 
on Athena (lower). The blue vertical dotted lines in Fig. [T] 
mark the orbital separations corresponding to periods of 
10®s, 10®s and lO^s for a 10 ®Mq MBH binary. From Fig. [U 
radial velocity shifts in the broad FeKa line due to MBH bi¬ 
naries with mass ratios q > 0.01 and separations of several 
hundred Vg could be measured and limits placed on the rate 
of hard MBH binary occurrence in the local Universe. Note 
that the secondary MBH is likely to have its own accretion 
disk. However, FeKa flux from the secondary in an strongly 
unequal mass binary is likely to be swamped by emission 
from the primary disk (luminosity oc MM), hence here we 
focus on emission from around the primary MBH. 


2.1 Precession of the secondary orbit 

As the MBH binary approaches merger, its orbit will tend to 
circularize due to emission of gravitational radiation, how¬ 
ever other factors may act against this effect (see below and 


^ Where the spin of the primary is important, Ant can be defined 
in terms of the precession of the orbital angular momentum of 
the secondary around the spin axis of the primary (AU) at fixed 
inclination and the precession of the periapse withi n the orbital 
plane (At^:), as AiU = Alo + cos(i)AO llMeri'itdl2013l ). 



Figure 1. The peak fractional change in energy SEJE from eqn. [2] 
as a function of mass ratio q and binary separation 05 . Hori¬ 
zontal dashed red lines denotes the planned energy resolution at 
6.40keV for Astro-H{'^ 6 eV) and Athena(2e\l). The blue vertical 
dotted lines mark the orbital separations correspond to oscilla¬ 
tion periods around the barycenter of 10 ®’®Fs respectively, for 
a binary of mass 10® Mq. For comparison, XMM-Newton EPIC 
has 5E/E ~ 0.025, making it insensitive to most MBH binary 
parameter space. 


e.g. ICuadra et al.l ll2009ll : iRoedig et al.l (l2012l l~). As the bi¬ 
nary orbit shrinks below at ~ 10®“^rg, general relativistic 
(GR) effects cause orbital changes that may be detectable 
in oscillations of the centroid energy of the broad FeKa line. 
Changes in uJ over time will have the effect of modulating the 
sinusoidal function of 5E/E in eqn. [21 The overall precession 
angle evolves as 

aj{t) = uj(0) -h [—CUdisk + OJlK + OJG-R + ojlt + ajQ]t (3) 

where aJ(0) is the initial angle value, Zudisk is the rate 
of (retrograde) precession of periapse due to the interac¬ 
tion between the secondary and the disk (analagous to the 
gravitational tugs between planets in the planetary case), 
lOlk is the rate of precession of periapse due to the inter¬ 
change of (e,i) due to the Lidov-Kozai (LK) mechanism, 
ZOgr is the (geodetic) rate of precession of periapse due to a 
Schwarzschild spacetime, lOlt is the precession rate of peri¬ 
apse due to Lense-Thirring (LT) frame-dragging by a spin¬ 
ning primary and lJq is the rate of precession of periapse 
due to the quadrupole moment of the non-spherical (spin¬ 
ning) primary. Note tha t LT frame-dragging effects can lead 
to inner disk warping (iBardeen fc PettersonI [19751 1. so the 
angle of the inner disk to the observers sightline can be in¬ 
dependent of the angle of orbital plane of the secondary {i 
in eqn 121. 

In planetary studies of the effect of planetary tugs on 
the radial velocity component of stars, aJ tends to evolve 
slowly with time and the effect of co-orbital planets (equiv¬ 
alent to our Zjdisk term above), dominates over other terms. 
However, in the case of a merging MBH binary, GR effects 
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Figure 2. SEjE over two periods of oscillation of the barycenter 
of a MBH binary inclined at 60° to the observers’ sightline (90° is 
edge-on). Red dashed lines represent the expected limits of SEJE 
for Astro-H. A 300M7aiooks exposure corresponds to one orbit 
at lOOrg around a Mj = Mi/IO^Mq SMBH. 


become important and we can approximate the precession 
of per lapse of th e secondary MBH orbit per 27r revolution 
as (lMerrittll2013l l 


Slu ~ G'kF 


1 ± 


4s 


g r’^/^cos(i) +jF(l-4cos^(i)) 


+5a;disk(4) 


where the term in square brackets is [5 cJgr + SUjut + ^oJq], 
F = TgllabD), s — ±[0,1] is the dimensionless spin pa¬ 
rameter of the primary MBH and (5Zjdisk = 2'Kqa),ll{'i -\-1), 
where is the mass ratio of distributed 

(disk) mass at r < 05 to mass Mi. The SlJgr = 6nF 
term dominates ^ZUdisk when D < {Svg/abqaD^^^, which 
holds if abjvg < Mi/M(j,^„j^) for modest e, or broadly 
if ab ^ lCd~'^rg for most realistic disk models. The con¬ 
tribution of spin terms (JZUspin = SUjtjt + SUjq) to the 
precession rate is small at wide binary separations, grow¬ 
ing from (5Zdspin/(5ZZJGR ~ 0.02 at 10®to 5aJspin/(5ZZJGR ~ 
0.2 at lOrg. We ignore the effect of the LK mecha¬ 
nism (JZUlk) since it requires a tertiary (disk or satellite) 
at high re lative inclination and more m assive than the 
secondary llFord. Kozinskv fc Rasiol l2000l l. However, since 
stalled SMBH are expected to occur o n ~pc scales in galactic 
nuclei (iMilosavlievic fc Merritlll200lll . it is possible that the 
LK mechanism could yield binary eccentricities e > 0 in ran¬ 
domly observed MBH binaries even at very small (~ lOOrg) 
separation. 

Fig.[2]shows SEjE for the centroid energy of the broad 
FeKa line as a function of a 2?! orbit over two periods of 
oscillation of the primary about the binary barycenter. Solid 
lines depict MBH binaries (with e = 0.01, 0.1) at ab = lOOrg 
and the black dashed line depicts an MBH binary at ab = 
200rg. A q — 0.01 binary at ab = lOOrg yields an oscillation 
that is not detectable with Astro-H (dashed red lines) but 
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Figure 3. As Fig. [2] except over 40 oscillation of the barycen¬ 
ter of an MBH binary. Curves with \SE/E\ ~ 0.02 show the 
envelope of oscillations for g = 1, = 500rg at eccentricities 

e = O.l(upper) and e = O.Ol(lower). Curves with \5E/E\ ~ 0.008 
show the envelope of oscillations (first three are depicted) for 
q = 0 . 1,05 = lOOrg at eccentricities e = 0.01,0.1. Curves with 
\&E/E\ ~ 0.001 show the envelope of oscillations (oscillations 20- 
23 are shown) for q = 0 . 01,05 = lOOrg at eccentricity e = 0.1. 
Red dashed lines represent the expected limits of SEJE for Astro- 
H. Precession of periapse (by ~ 11°/orbit at 05 = lOOrg) leads to 
periodicity in SEjE after ~ 33 orbits. A 1.2Ms exposure corre¬ 
sponds to 40 orbits at lOOrg around a Mq = MiM q SMBH. 


is detectable with Athena {5E/E ~ 3 x 10“"*). Even over 
only two periods of oscillation, the effects of GR precession 
on different binary configurations are starting to appear in 
the slight phase shift between the e = 0.01,0.1 curves and 
between the ab = 100, 200rg curves in Fig. [2l 

Fig. [ 3 ] shows an expanded version of Fig. [2] showing 
5E/E for the centroid energy of the broad FeKa line as a 
function of 27r orbit, over 40 orbits of the primary about the 
binary barycenter. The two periods of oscillation from Fig. [2] 
for g = 0.1, 05 = lOOr-g are shown in Fig.[3]as a part of a three 
period oscillation between ±0.008. The envelope for these 
oscillations assuming e = 0.01,0.1 is extended to 40 orbits, 
revealing the modulating effects of GR-dominated orbital 
precession of the secondary, showing up in oscillations of the 
primary around the barycenter. Also shown in Fig.^are the 
envelope of oscillations of the barycenter of a g = 0 . 01 , 05 = 
lOOrg, e = 0.1 MBH binary together with oscillations 20 — 23 
which will therefore not be detectable with Astro-H (red 
dashed lines), but will be detectable with Athena (assuming 
5E/E ~ 3 X 10“"^). Fig. [3] also shows the limiting case of 
the envelope of oscillations due to a g = 1 MBH binary at 
05 = SOOcg with e = 0.1, 0.01 (assuming only one MBH has 
an associated broad FeKa line for simplicity). Note that a 
1.2M6 Ms exposure corresponds to 40 orbits at lOOrg around 
a Me = Mi/10®Mq SMBH (see observing strategies below). 
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2.2 The broad and narrow FeKa line components 


Broad FeKo lines are ob served in several nearby active 
galac tic nuclei (AGN) le.g. iNandra et abll 19971 : iBraito et al.l 
I2OO7I ) and originate in fluorescent Fe deep in the gravi¬ 
tational pote ntial well of the central supermassive black 
hole (see e.g. Ife vnoMs fc Nowalj l200,ll . for a review). See 
iMcKernan et al' (l2013ll land references therein) for a dis¬ 
cussion of how the broad FeKa line component can be 
used to probe configurations of gas in the inner accre¬ 
tion disk and therefore test for possible MBH binaries. 
The narrow core of the FeKa li ne observed in AGN 
(e.g. lYaaoob fc PadmanabhanI 120041 ') originates in fluoresc¬ 
ing c old gas (Fe l-Fe xvil) far fro m the supermassive black 
hole llShu. Yaaoob fc Wandl201(]l ). We expect variations of 
the narrow core of the FeKa line to be uncorrelated on 
observing timescales with variations in the broad compo¬ 
nent and, more importantly, the narrow line centroid should 
not vary. Therefore, without knowing the absolute value of 
SE/E in the broad line in a given AGN, we can monitor its 
variation over time, by comparing the centroid of the narrow 
line with the oscillating centroid of the broad component. 

Fig. |4] shows a lOOks simulatec0 observation with the 
planned micro-calorimeter on Athena of an AGN at z = 0.01 
with a 2-lOkeV photon flux of 4.5 x 10“^^ergcm“^s“^. The 
simulated data shown (binned at ~ 12eV) correspond to 
a looks exposure with 5E/E equal to the average from 
the first positive oscillation for q = 0.1,= lOOrg in 
Fig. [3 (or SE/E = 4 X 10 The blue curve corresponds 
to the best-fit model to these data and the red curve shows 
the best-fit model to an equivalent lOOks observation with 
SE/E = —0.004, equivalent to the average SE/E for orbits 
0.5-1.0 in Fig. [21 The displacement between the two model 
fits by eye highlights how the oscillation of binaries around 
their barycenter will be easily detectable with Athena. The 
narrow FeKa line component (ignored in this model) should 
maintain the same centroid energy, allowing the relative shift 
in SE/E to be determined. The XSPEC toy model for the 
data in Fig.|3]was zphabs(zpowerlw-I-diskline) where the 
absorbing Hydrogen column was 10^^cm~^ and the power- 
law index was —1.8. The primary diskline for the simulated 
data had an inclination angle of 60° and emissivity index 
k = —2.5, as well as a line normalization of 0.1 with respect 
to the continuum to illustrate the effect of the radial-velocity 
shift. 


3 OBSERVING STRATEGIES 


From Fig. (2] a 27r barycenter oscillation period for a 
q = 0.1, Ub = lOOrg binary will take ~ lyear around 
a Ml = 10 ®Mq primary SMBH, but only ~ 8 hours 
around a Mi = 10® Mq primary SMBH. Since typical ex¬ 
posures of hundreds of ks (few days) are required to build 
up line fluxes with modern X-ray telescopes, techniques to 
search for MBH binaries in FeKa lines from AGN must 
depend on the moderately accurat e estimates of the cen¬ 
tral mass (e.g. Gulf^in et al. I I2OO9I : iMcKernan et alJI2OI0I : 
iKormendv fc hJ 2013h . Large mass (~ 10® M©) SMBH are 
the simplest to test for barycenter oscillations, since a MBH 


^ http://heasarc.gsfc.nasa.gov/cgi-bin/webspec 



Figure 4. A simulated lOOks observation with the Athena mi¬ 
crocalorimeter of an AGN at 2 = 0.01, with a 2-lOkeV flux of 
4.5 X 10~^^ergcm~ ^s~^, which is generated by a diskline model 
l|Fabian et al.lll989h (see text for details). Black data points cor¬ 
respond to simulated data where the average line centroid shift is 
SEjE = -1-0.004, which is the average of the first positive oscilla¬ 
tion in Fig. d and the blue solid line corresponds to the best-fit 
model to these data. The red solid line corresponds to the best 
model fit when SE/E = —0.004, or the average SE/E for the first 
negative oscillation in Fig. d (simulated data not shown for clar¬ 
ity). The 200ks oscillation period corresponds to a g = 0.1 binary 
at a(, = lOOrg with Mi = 6 X 10® Mq. 


binary will yield a monotonic change in the FeKa cen¬ 
troid energy between exposures over year (a^ ~ lOOrg) 
to multi-year (at > lOOrg) timescales. Moderate mass (~ 
10^ — 10® Mq) SMBH will yield an oscillation over a single 
long exposure- a SOOMrUiooks exposure corresponds to one 
orbit at ai, = lOOrg around a Mr = Mi/IO^Mq SMBH. In 
these cases, a division of the exposures into two or three seg¬ 
ments, at different phases, is sufficient to test for MBH bina¬ 
ries at moderate separations. The smallest mass (~ 10® Mq) 
SMBH will be the most difficult to test for MBH binaries, 
since a several hundred ks exposure might contain tens of 
oscillations. The oscillations will broaden the wings of the 
broad FeKa compon ent and models measur ing the spin of 
a single SMBH (e.g. iBrenneman et ahl 1201 if) will therefore 
over-estimate of the spin of the SMBH. To test for tight 
MBH binaries in low mass AGN, and constrain q, ab param¬ 
eter space, long exposures must be chopped into periodic 
odd-even sets on short timescales at different phases, that 
may be co-added to test for radial velocity oscillations and 
modulation by precession effects. Note that the effective area 
of Athena (which falls off at higher energies) may be well ex¬ 
ploited by AGN at modest redshifts. 
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4 CONCLUSIONS 

Analagous to the search for planets around stars using radial 
velocity shifts in stellar spectral lines, radial velocity shifts 
will be imprinted on the broad FeKa line due to oscillations 
of a massive black hole binary (MBHB) around its barycen- 
ter. Near future X-ray telescopes {Astro-H or Athena) will 
have the energy resolution to search for such oscillations in 
the broad FeKa line observed in active galactic nuclei. We 
demonstrate that oscillations of MBHB about their barycen- 
ter with mass ratios > 0.01, moderate binary separations 
(o-b ^ lOOOrg) and moderate eccentricities (e ~ 0.01 — 0.1) 
will be detectable in multiple exposures of the broad FeKa 
line over timescales appropriate for the MBH mass. The pe¬ 
riodic red- or blue-shift of the broad line energy centroid can 
be calibrated against the centroid of the narrow FeKa line. 
The magnitude of precession of the secondary MBH orbit is 
dominated by GR effects and we demonstrate that this will 
be detectable for modest mass ratios and separations. Fol¬ 
lowing Bayesian searches for radial vel ocity effects in stars 
due t o tugs by unobserved planets(e.g. iBrewer fc DonovanI 
bOliih . we can take a similar approach to finding periodic 
signal in integrated (and multi-epoch) exposures of broad 
FeKa lines in AGN. Observations of such oscillations will 
reveal the closest MBH binaries and brightest gravitational 
wave sources in the local Universe. 
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